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ABSTRACT 

Recent observational studies of intermediate-age star clusters (SCs) in the Large 
Magellanic Cloud (LMC) have reported that a significant number of these objects 
show double main-sequence turn-offs (DMSTOs) in their color-magnitude diagrams 
(CMDs). One plausible explanation for the origin of these DMSTOs is that the SCs 
are composed of two different stellar populations with age differences of ~ 300 Myr. 
Based on analytical methods and numerical simulations, we explore a new scenario 
in which SCs interact and merge with star-forming giant molecular clouds (GMCs) 
to form new composite SCs with two distinct component populations. In this new 
scenario, the possible age differences between the two different stellar populations 
responsible for the DMSTOs are due largely to secondary star formation within GMCs 
interacting and merging with already-existing SCs in the LMC disk. The total gas 
masses being converted into new stars (i.e., the second generation of stars) during 
GMC-SC interaction and merging can be comparable to or larger than the masses of 
the original SCs (i.e, the first generation of stars) in this scenario. Our simulations 
show that the spatial distributions of new stars in composite SCs formed from GMC- 
SC merging are more compact than those of stars initially in the SCs. We discuss 
both advantages and disadvantages of the new scenario in explaining fundamental 
properties of SCs with DMSTOs in the LMC and in the Small Magellanic Cloud 
(SMC). We also discuss the merits of various alternative scenarios for the origin of the 
DMSTOs. 
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1 INTRODUCTION 

Recent photometric studies of intermediate-age SCs in the 
LMC have discovered objects exhibiting extremely unusual 
main-sequence turn-offs (MSTOs) in their CMDs. The first 
indications that such clusters might exist date back several 
years to observations obtained with terrestrial facilities - 
Bertelli et al. (2003) demonstrated that the LMC cluster 
NGC 2173 apparently possesses an unusually large spread in 
colour about its MSTO, while Baume et al. (2007) obtained 
a similar result for NGC 2154. 

More recent studies based on deep precision photometry 
from the Hubble Space Telescope (HST) Advanced Camera 
for Surveys (ACS) have revealed the truly peculiar nature of 
many intermediate-age LMC clusters. Most striking are the 
relatively massive clusters NGC 1846, 1806, and 1783 which 
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possess double MSTOs (hereafter referred to as DMSTOs for 
convenience) on their CMDs (Mackey & Broby Nielsen 2007, 
MBN07; Mackey et al. 2008; Milone et al. 2008; Goudfrooij 
et al. 2008). Despite this, the remaining CMD sequences 
for these clusters (i.e., their red-giant branches, subgiant 
branches, main sequences, and red clumps) are all extremely 
narrow and well-defined. This suggests that none of the clus- 
ters possesses a significant line-of-sight depth or internal dis- 
persion in [Fe/H], or suffers from significant differential ex- 
tinction. The simplest explanation is that the observed DM- 
STOs in each of these clusters represent two distinct stellar 
populations with differences in age of ~ 200 — 300 Myr. It 
should be noted, however, that apart from the study of Muc- 
ciarelli et al. (2008) who found no significant star-to-star dis- 
persion in [a/Fe] for 6 stars in NGC 1783, no constraints on 
the possibility of internal variations in chemical abundances 
other than [Fe/H] presently exist for these SCs. 

Milone et al. (2008) studied a large sample of 16 
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intermediate-age LMC clusters, and found that 11 of these 
possess CMDs exhibiting MSTOs which are not consistent 
with being simple, single stellar populations. Of these 11, 
four clearly show DMSTOs (the three described above plus 
NGC 1751) while the remainder possess more sparsely popu- 
lated CMDs that make the precise morphologies of their pe- 
culiar MSTOs difficult to ascertain. All are consistent with 
being DMSTOs, or alternatively with possibly being more 
smoothly distributed intrinsic broadenings. Nonetheless, all 
11 clusters again exhibit very narrow sequences across the 
remainder of their CMDs, suggesting that their peculiar 
MSTOs are due to internal age dispersions of ~ 200 — 300 
Myr. Milone et al. (2008) measured the ratio between the up- 
per and the lower MSTO population in NGC1846, NGC1806 
and in the less-populated NGC1751 and found that the in- 
ferred young populations may comprise up to ~ 70 per cent 
of the stars in the central regions of these SCs. 

In addition to the above objects, several younger LMC 
clusters have been identified as possibly possessing irregular 
CMDs (e.g., Santiago et al. 2002, Gouliermis et al. 2006). 
It is further worth noting that one intermediate-age SMC 
cluster, NGC 419, also likely possesses a DMSTO (Glatt et 
al. 2008). 

The results described above suggest that, contrary to 
expectation, DMSTOs may be a common feature among 
intermediate-age SCs of the Magellanic Clouds (MCs). If 
the origin of such DMSTOs is indeed closely associated with 
an age difference between a first and second generation of 
stars in these SCs, the above observations raise the following 
three key questions: (i) why the second generation of stars in 
a SC with a DMSTO can be formed about 300 Myr after the 
formation of the first generation of stars, (ii) why the total 
mass of the second generation can be comparable or even 
larger than that of the first, and (iii) how the two different 
populations can now be in the same SC. Although the origin 
of the observed multiple stellar populations in the Galactic 
globular clusters (GCs) has been observationally and theo- 
retically discussed in terms of GC formation scenarios (e.g., 
Piotto 2008), the above fundamental questions so far have 
not been discussed by theoretical models of SC formation in 
the LMC. 

The purpose of this paper is to explore a new scenario 
in terms of the above three fundamental questions on the 
origin of the observed DMSTOs. In the new scenario, an 
already formed SC (i.e., the first generation of stars) can 
merge or interact with a CMC so that a second generation 
of stars is formed in the CMC. The second generation can 
then merge and mix dynamically with the SC to form a new 
SC with two stellar populations of different ages. Thus, the 
origin of the DMSTOs in SCs of the LMC results from merg- 
ing between SCs and GMCs in the new "GMC-SC-merger" 
scenario. The typical age difference (i gap ) of the two popu- 
lations in SCs required for explaining the DMSTOs corre- 
sponds to the typical time scale of merging between SCs and 
GMCs of the LMC in the new scenario. We mainly discuss 
(i) whether star formation is likely to be enhanced during 
GMC-SC merging, (ii) the mass fraction of the second gener- 
ation of stars relative to the first in a SC and its dependence 
on parameters of GMC-SC merging, and (iii) time scales of 
GMC-SC merging in the LMC. 

The plan of the paper is as follows: In the next section, 
we describe our numerical models for GMC-SC merging in 



the LMC. In §3, we present numerical results mainly con- 
cerning the final distributions of gas and new stars in the 
remnants of GMC-SC mergers. In §4, we discuss the origin 
of the DMSTOs of SCs in a more general way. In this sec- 
tion, we also discuss the advantages and disadvantages of 
two alternative scenarios in explaining the physical proper- 
ties of SCs with DMSTOs. We summarize our conclusions 
in §5. 



2 THE MODEL 

The most important question in the present study is whether 
the masses of the second generations of stars (i.e., new stars) 
in the remnants of GMC-SC mergers can be as large as those 
of the first (i.e., stars initially in SCs). For the masses of the 
second generations to be comparable to those of the first, a 
large amount of gas initially in GMCs needs to be accumu- 
lated within merging SCs for star formation. Although we 
do investigate "star formation processes" within GMCs, the 
adopted model for star formation is rather idealized and phe- 
nomenological, and does not allow us to investigate subpc- 
scale real star formation processes (e.g., the formation of 
individual stars within the cores of small molecular clouds). 
We therefore consider that here it is more important for us 
to (i) investigate how much gas in GMCs can be captured 
in SCs during GMC-SC merging for various different model 
parameters for the merging and (ii) suggest that the merger 
remnants (i.e., new SCs) with large fractions (~ 0.3 — 0.5) of 
gas can finally become SCs with two different populations 
with comparable masses. 

We adopt the same numerical code as adopted in our 
previous numerical simulations on the hydrodynamical evo- 
lution of interstellar gas in interacting and merging galax- 
ies (Bekki et al. 2002) and the formation of star clusters 
from GMCs in galaxies (Bekki & Couch 2003). Although we 
can investigate gas dynamics in merging GMCs, the adopted 
TREESPH code does not allow us to precisely describe long- 
term stellar dynamics such as two-body relaxation processes 
in SCs. In future work we will discuss the long-term evolu- 
tion of the remnants of GMC-SC mergers by using appro- 
priate numerical codes such as NB0DY4 adopted in our pre- 
vious simulations of SC formation (Hurley & Bekki 2008, 
HB08). In order to understand more clearly the hydrody- 
namical evolution of GMC-SC mergers, we do not include 
the external tidal field of the MCs in the present study. 

A SC in a GMC-SC merger is represented by a Plum- 
mer model with mass and size represented by M B and R B , 
respectively. The scale length of the Plummer model is fixed 
at 0.2i? s for all our simulations. A GMC is assumed to be a 
Bonncr-Ebert isothermal sphere (Ebert 1955; Bonner 1956) 
with a mass of M s , a size of i? g , and a temperature of T g . 
The radial distribution of a Bonner-Ebert isothermal sphere 
is determined by solving the differential equations for a given 
T g of the sphere. Guided by the observed relation between 
mass densities and sizes of GMCs discovered by Larson's 
(1981) relation and the observed typical mass and size of 
GMCs in the Galaxy (e.g., Solomon et al. 1979), we use the 
following R s — M g relation; 
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Table 1. Model parameters for numerical simulations of GMC-SC mergers 

Model no a M s " R s c y g " |ug| e Star formation f Comments 

Ml 5.0 10 6.0 2.2 NO the standard model 

M2 5.0 10 6.0 2.2 YES the star formation model 

M3-32 (1.0, 20.0) (5, 25) (0, 23.8) (0, 12.8) NO & YES 



a For models M3 - M32, the ranges of model parameters investigated in the present study are shown in parentheses: the left and the 
right numbers are the minimum and maximum values. 

The total mass of stars in a SC in units of 10 4 Mq . The total gas mass m a CMC is 5 X 1O 4 M for all models. 
c The size of a SC in units of pc. The size of a CMC is set to be 11.9 pc for all models in the present study. 

" Initial y-position of a CMC with respect to the center of a SC. This corresponds to the impact parameter of the GMC-SC collision. 
e The absolute magnitude of the initial velocity of a CMC along the rr-axis in units of km s — 1 , with respect to the center of a SC. For 
all models, « g $C 0. 

f "YES" ("NO") means that the star formation model is (is not) included in the simulation. 



We investigate models with M g = 5 x 1O 4 M , R g = 11.9 
pc, and T g =2.6 K in the present study. Since the isother- 
mal Bonner-Evert sphere is adopted, one of key parameters 
determining the gas dynamics of GMC-SC merging is the 
mass-ratio (s g = M s /M s ) of a CMC to a SC in a given 
model: if we derive parameter dependences of numerical re- 
sults on s g for the above M g , these dependences can be 
applied for other M s . 

The initial locations of the SC and the CMC in the 
GMC-SC merger are set to be (a; B1 y s , z B ) and (x g , y g , z g ), 
respectively. For convenience, (x s , y s , z s ) =(0,0,0) so that the 
initial position of the GMC with respect to the SC can be 
more clearly understood. The initial velocities of the SC and 
the GMC in the x-, y-, and z-directions are set to be (u s , v s , 
w s ) and (lt g , tig, W g ), respectively, and (u s , v s , W s ) = (0,0,0) 
for simplicity. The orbital plane of the GMC-SC merger is 
coincident with the x-y plane for all models in the present 
study (i.e., z s = z g = w g = 0). We adopt x g = 2R g for 
all models so that GMC and SCs are not initially in direct 
contact. The initial direction of the velocity vector of the 
GMC is only a parameter for the orbit of the GMC owing to 
the spherically symmetric distributions of the model GMCs 
and SCs. We therefore assume that the initial direction is 
parallel to the i-axis toward the negative x (i.e., u g < and 
v g = 0). Thus j/ g and u g are free parameters that determine 
the orbit of the GMC with respect to the SC. 

To summarize, the key model parameters in the present 
simulations are M B (or s g ), R s , y g (i.e., the impact parame- 
ter) , and u g . 

For all models the total number of particles (TV) is 
40 000, half of which is for stars, and the initial gravita- 
tional softening length is 0.24 pc. We choose these particle 
numbers in order to facilitate investigation of the long-term 
evolution of merger remnants in the external LMC potential 
in our future numerical studies using NB0DY4, as adopted 
in our previous simulations (HB08): models of SCs using 
these codes can be undertaken within a reasonable time scale 
for N ~ 10 4 - 10 5 without excessive numerical costs (e.g., 
D'Ercole et al. 2008; HB08). 

Owing to the limited size and mass resolution of the 
present simulations, subpc-scale physical processes of star 
formation cannot be investigated in the present study. We 
do, however, try to investigate whether the gas accumulated 
in the central regions of SCs after GMC-SC merging can be 
sufficient to form new stars by adopting a simple prescrip- 



tion for star formation. In the models with "star formation" , 
a gas particle is converted into a collisionless new stellar 
one if the gas particle meets the following conditions: (i) 
the dynamical time scale of the SPH gas particle is shorter 
than the sound crossing time , and (ii) the gas is converging 
(i.e., Vv < 0, where v is the velocity vector of the gas par- 
ticle). These two conditions mimic the Jeans gravitational 
instability for gaseous collapse. The adopted model for star 
formation is rather idealized so that "new stars" need to 
be interpreted as possible formation sites of stars in GMCs. 
In the rest of this paper, the stars initially in SCs are just 
referred to as "stars" to distinguish them from new stars 
formed from gas. 

Although we investigate many models spanning differ- 
ent combinations of the four key parameters, we mainly 
present the results of the "standard model" which shows 
the typical behavior of gas accumulation in SCs merging 
with GMCs. We also briefly describe the results of the "star 
formation model" in which the simple treatment of star for- 
mation is included and which as parameter values exactly 
the same as those adopted in the standard model. 

The parameter values of the standard model and the 
ranges of parameters investigated for an additional 31 mod- 
els are summarized in Table 1. The standard model and the 
star formation model are represented as Ml and M2, re- 
spectively. Parameter values in models other than Ml are 
not described in the table in order to avoid consuming an 
unnecessarily large space for the description. In the present 
study we mainly show the dependences of the gas mass frac- 
tions within a radius of 10 pc (/ g ) in merger remnants on 
y g , u g , M s , and R a . 

We note that SC-GMC merging does not occur in some 
of our models with large y s and |w g |; however, gas-transfer 
between SCs and GMCs is still possible in these cases. We 
also describe the results of these models briefly in the present 
paper. It should be stressed that the present study is the first 
step toward understanding the possibly complicated interac- 
tions between SCs and GMCs in the MCs: the present mod- 
els are rather idealized in some respects, in order to grasp 
the essential ingredients of GMC-SC merging processes. We 
plan to investigate a more fully-consistent and sophisticated 
model for GMC-SC merging in our future studies. 
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Figure 1. Time evolution of stellar (magenta) and gaseous 
(cyan) distributions projected onto the x-y plane (upper four pan- 
els) and the x-z plane (lower four panels) for the standard model 
Ml. The time T shown in the upper left corner of each panel is 
given in units of Myr. Note that a gaseous disk can be seen in 
the inner region of the SC, and the center of the disk seems to 
deviate from the center of the SC. 



3 RESULTS 

3.1 The standard model 

Fig. 1 shows how gas in a GMC evolves during major merg- 
ing between the GMC and a SC with s g = 1 in the standard 
model Ml. During the off-center collision between the GMC 
and the SC, the central part of the GMC is tidally com- 
pressed to form a compact gaseous core. Owing to dynami- 
cal friction during merging, the core can sink into the central 
region of the SC to finally form a flattened gaseous spheroid 
(or gas disk) within the SC. The SC is not destroyed by this 
merging and thus retains its initial spherical shape, though 
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Figure 2. The time evolution of stellar (solid) and gaseous (dot- 
ted) masses within the central lOpc of the SC in the standard 
model Ml. 
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Figure 3. The cumulative stellar (solid) and gaseous (dotted) 
masses within radius R from the center of the SC at T = 11.0 
Myr in the standard model Ml. 

it loses about 30% of its initial mass owing to tidal stripping. 
Hence, the newly formed SC has a very compact, flattened 
gaseous disk in its central region: further star formation from 
this central gas can result in two distinct stellar populations 
in the SC. 

About 30% of the gas in the GMC is rapidly stripped 
during merging to form a very diffuse gaseous halo around 
the new SC. This stripped gas cannot be accreted onto the 




Figure 4. Dependences of / g (i.e., gas mass fraction within 
a radius of 10 pc in merger remnants) on y g (upper left), |u g 
(upper right), M B (lower left), and R B (lower right). The dotted 
and dashed lines represent s g = 1 and 2.5, respectively. 
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Figure 5. Final distributions of stars (magenta), gas (cyan), 
and new stars (yellow) projected onto the x-y plane (left) and 
x-z plane (right) for the star formation model M2, in which the 
model parameters are exactly the same as those adopted in the 
standard model, but star formation is included. 

SC after merging even in this isolated model, and thus is 
likely to be returned back to interstellar medium (ISM) of 
the LMC if interaction between the stripped GMC gas and 
the ISM in the LMC is included in future simulations. Fig. 
2 shows that the gas mass fraction within the central 10 pc 
of the new SC (/ g ) increases as the SC loses its stellar mass 
and gas is accreted into the SC. Fig. 2 also shows that / g 
can be finally as large as 0.5 at T — 11 Myr. This implies 
that if all of the gas can be converted into new stars, the 
mass ratio of the second generation of stars to the first can 
be as large as 1. We discuss later important implications of 
this result in terms of the origin of the DMSTOs. 

Fig. 3 shows that the final cumulative radial mass dis- 
tributions (M(< R), where R is the distance from the mass 
center of the newly- formed SC) of gas and stars at T = 11 
Myr are quite different in the sense that the gas shows a 
more compact distribution than do the stars. The half-mass 
radius is 3.3 pc for the stars and 1.0 pc for the gas, which re- 
flects that the collapsed central part of the GMC is directly 
transfered to the center of the SC owing to dynamical fric- 
tion. Therefore the gas mass fraction is significantly larger in 
the inner region of the new SC, which means that this region 
of the new SC would be dominated by young stars if star 
formation were to subsequently occur in the accumulated 
gas. 



3.2 Parameter dependences 

The dependences of the final properties of SCs on a parame- 
ter are investigated by changing the values of the parameter 
and by fixing those of other parameters. Dependences of the 
present numerical results on the four key parameters (j/ g , 
|ti g |, M s , and R s ) are described as follows: 

(i) As shown in Fig. 4, f s does not depend strongly on 
y g (i.e., impact parameter) as long as GMCs merge with 
SCs (i.e., t/g < 12 pc). However / g can be very small (< 0.1) 
for models with large y g , in which no merging occurs and 
only small fractions of gas are transferred to SCs during 
the GMC-SC tidal interaction. For example, the model with 
y g = 18 pc and other parameter values being the same as 
those in the standard model shows f s — 0.07, corresponding 
to a gas mass within the central 10 pc of 3.2 x 1O 3 M0. 

(ii) The models with high |u g | (> 7 km s _1 ) show very 
small / g (< 0.1), because high-speed encounters between 



GMCs and SCs can again prevent GMCs from merging with 
the SCs (see Fig. 4). There is no strong dependence of f s 
on |« g | for models with |u g | ^ 7 km s _1 in which GMC-SC 
merging does occur. Models with |u g | > 12 km s _1 show 
very small / g (< 0.01) even if y g = pc. 

(iii) GMCs can either destroy SCs or tidally strip most 
of their stellar envelopes, if the masses of the GMCs are 
significantly larger than those of the SCs interacting with 
the GMCs. Therefore the compact gas cores of GMCs cannot 
sink into the central regions of SCs in models with large s g 
(= M S /M B ). As a result of this, the models with large s g (> 
2.5) show small / g in the merger remnants for a reasonable 
set of other model parameters (see Fig. 4). It should be, 
however, stressed that the models with y g — 0, small |u g | 
(< 4 km s _1 ), and large s g (> 2.5) show large / g (> 0.7), 
because in these situations the SCs become trapped in the 
central cores of the GMCs. It might not be appropriate to 
call these very gas-dominated merger remnants "SCs". 

(iv) In models with smaller s g (larger M s ), SCs can 
accumulate a larger amount of gas from merging GMCs. 
However, in these models it is not possible for s g to be par- 
ticularly large, owing to the initially larger stellar masses in 
these models. Therefore, there is an optimum s g for which 
/ g is a maximum for a given set of model parameters. As 
shown in Fig. 4, the major merger model (i.e., s g = 1) shows 
the largest / g , which implies that GMC-SC mergers with 
comparable masses can form larger fractions of the second 
generations of stars after merging. 

(v) The models with larger _R S show larger / g for R s ^ 
20 pc. This is mainly because the final stellar masses within 
the central 10 pc for SCs with larger R s (and thus lower 
stellar densities) tend to be small owing to the more effi- 
cient stripping of stars during GMC-SC merging: the actual 
masses of gas within the merger remnants do not depend 
strongly on Rs for R s ^ 20 pc. The models with large R s 
(e.g., R s = 25 pc) show very small / g owing to the total 
destruction of SCs during GMC-SC merging. 

3.3 Star formation 

Owing to our adopted prescription for star formation in 
GMCs, new stars are formed preferentially in rather high- 
density regions of GMCs. New stars can therefore be formed 
during the accumulation of compact and flattened gas 
spheres in the central regions of GMCs due to strong tidal 
compression of the GMCs by SCs. The total mass of new 
stars within a merger remnant (i.e., newly-formed SC) de- 
pends on how much gas is accumulated into the SC during 
merging. Therefore the parameter dependences of mass frac- 
tions of new stars in merger remnants for models including 
star formation are very similar to those of / g for models 
without star formation. Thus we here describe only the re- 
sults of the standard star formation model M2. 

Fig. 5 shows that the new stars in the merger remnant 
show a compact, disky distribution and thus have quite dif- 
ferent spatial distributions from the original stars: there ap- 
pears to be a SC within a SC. This "SC-within-SC" ap- 
pearance is one of key characteristics of the remnants of 
SC-GMC mergers in the present study. The total masses of 
stars, gas, and new stars are 3.7 x 10 4 Mq, 1.1 x 10 4 Mq, and 
2.0 x 10 4 Mq, respectively, which means that the mass frac- 
tion of new stars among all stars is 0.35. The half-mass radii 
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for stars, gas, and new stars are 4.0 pc, 3.8 pc, and 1.2 pc, 
respectively, which clearly demonstrates that the new stars 
have a much more compact distribution than do the origi- 
nal stars. The simulated flattened structure of new stars is 
rather similar to the inner flattened "core" observed in M15 
(e.g., van den Bosch et al. 2006). 

Fig. 6 shows that the star formation rate (SFR) becomes 
as high as 0.05 M yr -1 in the star formation model M2 
when the compact gas sphere of the GM C is sinking into the 
SC owing to dynamical friction. The SFR is quite peaked 
so that there is little age dispersion among the new stars 
formed during GMC-SC merging. The model with larger y s 
(=11.9 pc) and thus a weaker tidal field from the SC shows 
a smaller peak SFR (~ 0.02Mq yr -1 ), in the later phase 
of GMC-SC merging. These results clearly suggest that the 
strength of the tidal field of a SC merging with a CMC is 
a key factor that can determine the peak SFR and the star 
formation efficiency of the GMC-SC merger. Fig. 7 confirms 
that the inner regions of the merger remnant are dominated 
by new stars, as suggested in Fig. 5. 

The derived nested spatial distributions (i.e., SC- 
within-SC appearance) of the simulated SCs are apparently 
not observed in the intermediate-age SCs of the LMC which 
have been tested (e.g., MBN07). Furthermore, observations 
have not yet discovered inner disky structures of SCs in 
the LMC, though the SCs themselves have rather flattened 
shapes in general (e.g., van den Bergh 2000). To resolve these 
inconsistencies in structural properties between the simu- 
lated and observed SCs requires that long-term dynamical 
relaxation processes play a role in dynamically mixing old 
and new components and thus in forming the canonical ra- 
dial density profiles observed in the present SCs (e.g., King- 
type profiles). 

The median relaxation time in a SC like NGC 1846 
should be ~ 2 Gyr (see for example the models in Mackey et 
al. 2008b), and a factor of a few shorter than this in the very 
centre. Therefore, there has conceiveably been enough time 
for dynamical relaxation processes to influence and modify 
cluster structure. It should also be noted that the second 
generation of stars in a composite SC formed as described 
above will initiate a new phase of violent relaxation in the 
SC as the most massive members evolve quickly and die, and 
any residual gas is expelled. This can be very good at dy- 
namically mixing the cluster on a relatively short time scale 
(see e.g., Meylan & Heggie 1997 and references therein). It 
is likely that these evolutionary processes might alter signif- 
icantly the centrally concentrated spatial distribution of the 
second generation of stars, as well as the observable number 
fraction of the first and second generations of stars in the 
centres of new SCs. 



4 DISCUSSION 

4.1 Advantages and disadvantages of the scenario 

The observed apparently clear distinction between bluer and 
redder MSTOs in some SCs of the LMC (e.g., MBN07; 
Mackey et al. 2008; Milone et al. 2008) implies that there 
were two bursts of star formation rather than prolonged pe- 
riods of star formation in these SCs. One advantage of the 
present scenario is therefore that it can naturally explain 
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Figure 6. The time evolution of star formation rates (SFRs) 
in units of yr for star formation models with j/ g = 6.0 pc 
(solid) and y g = 11.9 pc (dotted). For these models, gas can be 
converted into new stars if the gas particles satisfy the requisite 
conditions for star formation as described in the main text. Pa- 
rameter values other than y g for these two models are exactly the 
same as those of the standard model Ml. 
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Figure 7. The same as Fig. 3 but for stars (solid) and new stars 
(dotted) in the standard star formation model M2. 

the possible two burst epochs of star formation: one is dur- 
ing the original formation of SCs and the other is during 
GMC-SC merging. Furthermore, the new stars formed dur- 
ing the GMC-SC merging are clearly bound to the original 
SC, which naturally explains how two different stellar popu- 
lations may coexist in some SCs. The scenario does not need 
to consider a physical mechanism by which interstellar gas 
can be accumulated within SCs and then form stars within 
the SCs. 

The present scenario predicts significant age differences 
between the first and the second generations of stars in SCs 
owing to the late formation of stars triggered by GMC- 
SC merging/interaction. A key question here is therefore 
whether the age difference can be similar to ~ 3 x 10 yr 
(igap), which may be required to explain the origin of the 
DMSTOs (e.g., MBN07; Mackey et al. 2008; Milone et al. 
2008). The time scale of a GMC-SC merger event (t m ) in 
the LMC can be estimated as follows (e.g., Makino & Hut 
1997): 



where n g , a, and v are the mean number density of the 
GMCs, the geometrical cross section of a CMC, and the 
relative velocity between a GMC and a SC. If we consider 
that the LMC is a uniform disk with a disk radius of r<j (~ 5 
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kpc), a scale height of z d , and the number of GMCs being 
equal to N g for simplicity, then 

- = (3) 
We assume here that SCs are more compact than GMCs 
with comparable masses and thus that a in the i m estimation 
should be the cross section for GMCs. If we consider that 
GMCs have spherical shapes, then we can use the following: 

a = Tvr g 2 , (4) 

where r g is the size of a GMC (i.e., r g = 0.5R g ). 

We are interested in the SCs of the LMC a few Gyr ago 
when the SCs with DMSTOs were formed. The LMC would 
have a larger gas mass fraction and thus a larger number 
of GMCs in the past. Recent NANTEN observations have 
found 272 GMCs with masses larger than 2 x 10 4 M Q in the 
LMC (Fukui et al. 2008) ; however the number of GMCs a few 
Gyr ago may be significantly larger than that in the present 
LMC. If we consider (i) that v is similar to the observed 
velocity dispersion of stars in the LMC (e.g., van der Marel 
2002), (ii) that the past iV g is larger by a factor of 2 than 
that of the present LMC, and (iii) that z d is about 10% of 
the observed radial scale length of the LMC (=1.5 kpc; van 
den Bergh 2000), then t m can be estimated as follows: 



^-^Op ( 150pc )( 15pc ; 



- 20kms 



~) Gyr. (5) 



It should be stressed here that the above i m is for the average 
number density of GMCs with masses larger than 2 x 1O 4 M0 
over the entire LMC region: £ m can be much shorter than 
the above in some regions with locally high number densities. 
The above value of i m is much larger than i gap (~ 0.3 Gyr) 
required for explaining the DMSTOs. However, i m can be 
as small as i gap in the LMC if we adopt a higher N g , an r g 
larger than equation (1) describes, and a thinner gas disk as 
follows: 
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It remains unclear which of the above two estimations is 
more reasonable and realistic in the LMC a few Gyr ago. In 
addition, owing to differences in the structural and kinemat- 
ical properties of gas and stars between the LMC and the 
SMC (i.e., N g , Zd, r g , ra, and v), tm might well be different 
between the MCs. 

A SC can merge with a GMC located within a distance 
d g from the SC within the following time scale: 

r ) " G yr- (7) 
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This means that t m can be similar to t gap for some local 
regions with surface number densities (E g ) of GMCs similar 
to 73 kpc -2 . Such high S g would be possible in strong tidal 
arms and in the inner region of the LMC. It is also possible 
that a group of GMCs (or GMC associations) could have 
such a high E g in the LMC. Thus, £ m can be as small as 
~3x 10 s yr (thus ~ i gap ), although it is possible that i m 
varies significantly between different regions of the LMC a 
few Gyr ago. 

It appears that the observed possible age differences 
between two populations is an order of 10 s yr: no SCs with 
DMSTOs appear to show implied age differences of an order 



of 10 9 yr (see Santiago et al. 2002 for a possibly exceptional 
case of NGC 1868). Therefore one possible problem in the 
present scenario is that the time scales of GMC-SC merging 
can in principle be significantly longer than 10 8 yr for some 
SCs: the scenario needs to explain why there are apparently 
no intermediate-age SCs with implied age differences be- 
tween the two component stellar populations being as large 
as 1 Gyr. A possible solution of this potential problem is if 
the time scales of GMC-SC merging in the LMC grew longer 
in the more recent past - that is, as the LMC evolved with 
time GMC-SC merging became less likely because the num- 
ber density of GMCs and and the typical GMC size became 
lower and smaller, respectively, owing to rapid gas consump- 
tion by ongoing star formation in the LMC. Thus it would 
be possible that intermediate-age SCs cannot merge with 
GMCs much later (~ 1 Gyr) than their formation in the 
LMC. 

Furthermore, it would be equally possible that the time 
scales of GMC-SC merging can be significantly shorter than 
10 yr owing to locally high number densities of GMCs in 
the vicinity of some SCs. Although SCs which experience 
merging with GMCs much less than 10 8 yr after their for- 
mation can still have two distinct component populations, 
they would not clearly show DMSTOs in their CMDs due to 
the much smaller age differences between them. Therefore 
the possible presence of SCs with age differences significantly 
smaller than 10 8 yr would be still consistent with observa- 
tions. Merging of SCs with GMCs less than 10 8 yr after 
their formation might well be responsible for the observed 
binary/multiple SCs with small age differences (e.g., NGC 
1850, Leon et al. 1999). 

GMC-SC merging is highly likely to occur between SCs 
and their local neighborhood GMCs in the LMC: cloud com- 
plexes (or "superclouds" ) with the typical mass of 1O 7 M 
(Elmegreen 1987) would be the progenitors for the GMC- 
SC mergers. Therefore differences in chemical abundances 
between merging SCs and GMCs should be quite small un- 
less there exist significant radial and azimuthal abundance 
gradients in the LMC. Owing to the time lag of ~ 3 x 10 s yr 
between SC formation and SC-GMC merging in the GMC- 
SC merger scenario, GMCs can be chemically polluted by 
AGB stars in SCs merging with the GMCs. Therefore, it is 
quite possible that the second generations of stars formed 
from GMCs could have different abundances in light ele- 
ments from those in the first generations initially in the 
SCs. Ongoing spectroscopic observations on the abundance 
properties of SCs with DMSTOs will soon reveal whether or 
not abundance inhomogeneities exist in these objects. We 
thus plan to investigate the chemical evolution of GMC-SC 
mergers in order to compare between the predicted and the 
observed degrees of inhomogeneity in light element abun- 
dances for SCs with DMSTOs in the LMC. 



4.2 Alternative scenarios 

4-2.1 Self-pollution by AGB stars 

It would be possible that stellar ejecta from AGB stars in 
a SC can be accumulated in the central region of the SC 
so that the ejecta can be used for the formation of the sec- 
ond generation of stars. In this "self-pollution" scenario, the 
second generation of stars needs to form well after the re- 
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moval of stellar ejecta from Type II supernovae that can 
cause a significant spread in heavy element abundances. Re- 
cent numerical simulations of GC formation based on this 
self-pollution scenario have shown (i) that star formation for 
the second generation can start soon after ejection of stellar 
winds from massive AGB stars, and (ii) that star formation 
can continue gradually within GCs over a period of ~ 100 
Myr (D'Ercole et al. 2008). 

Thus the self-pollution scenario has serious difficulties 
in explaining (i) why the time difference between the for- 
mation epochs of the first and second generations of stars 
is typically ~ 300 Myr, and (ii) why most stars in the sec- 
ond generation formed almost simultaneously. It is likely 
that the numerical simulations by D'Ercole et al. (2008) do 
not describe so precisely the star formation histories within 
SCs owing to the adopted rather idealized models for hydro- 
dynamics and star formation within SCs. Therefore, future 
theoretical studies based on more sophisticated numerical 
simulations are certainly worthwhile to confirm whether the 
self-pollution scenario has really the above serious problems. 

4-2.2 Star cluster merging 

It is a well known observational fact that a significant frac- 
tion of SCs in the LMC are binary or multiple clusters (e.g., 
Bhatia & Hatzidimitriou 1988; Bhatia et al. 1991; Dieball et 
al. 2002). This fact implies that merging of binary SCs can 
form single SCs with distinct stellar populations, if there are 
initial differences in age between the two original clusters. 
This merger scenario needs to explain why age differences 
between two merging clusters can be similar to 3 x 10 s yr for 
some binary SCs. Observational studies on possible age dif- 
ferences in binary SCs in the LMC have demonstrated that 
the components of binary or multiple SC systems typically 
appear to be small, hence implying that all SCs in a given 
bound system are generally formed simultaneously or over 
a very short time scale (e.g., Dieball et al. 2002). Therefore, 
the merger scenario appears to have a serious problem in ex- 
plaining the frequently occurring age difference of ~ 3 x 10 8 
yr observed for the DMSTOs in intermediate-age MC clus- 
ters. 

It should, however, be stressed that there are some bi- 
nary SCs in the LMC with age differences up to ~ 0.5 Gyr 
(e.g., SL356-357 pair, Leon et al. 1999): it would be possible 
that merging of the two component SCs in such systems can 
happen in the LMC. If this is the case, the merger scenario 
needs to explain (i) why SC merging apparently happens 
preferentially for those systems which have age differences 
of ~ 300 Myr and (ii) why the mass-ratios of the younger 
SCs to the older ones are comparable to and larger than 1. 
It is currently unclear whether such preferential SC merging 
can happen in the history of the LMC. 

4.3 DMSTOs only for intermediate-age SCs? 

The age distribution of the LMC SCs shows a gap extend- 
ing from 3 to 13 Gyr - with only one cluster in this age 
range - suggesting that a second epoch of cluster forma- 
tion started abruptly in the LMC about 3 Gyr ago (e.g., Da 
Costa 1991; Geisler et al. 1997; Rich et al. 2001; Piatti et al. 
2002). SCs with DMSTOs are observed to be - 1.5-2.5 Gyr 



old (Mackey et al. 2008; Milone et al. 2008) and thus were 
formed just after the "age-gap" period. So far it remains ob- 
servationally unclear whether the oldest GCs in the LMC, 
which formed before the commencement of the age gap, also 
show DMSTOs. This raises the following two questions: (i) 
whether the origin of the SCs with DMSTOs is closely asso- 
ciated with some specific formation processes at the reacti- 
vation of SC formation a few Gyr ago, and (ii) whether SCs 
can have DMSTOs irrespective of their formation epochs 
owing to a general physical mechanism responsible for the 
DMSTO formation. 

Recent numerical simulations have shown that tidal in- 
teraction between the LMC and the SMC can dramatically 
increase (by a factor of ten) the cloud-cloud collision rate 
in the LMC during the interaction (e.g., see Fig. 1 in Bekki 
ct al. 2004a). This result implies that if a high cloud-cloud 
collision rate means a high GMC-SC collision/merging rate 
in the interacting MC system, the origin of SCs with DM- 
STOs might be closely associated with the commencement 
of strong tidal interaction between the MCs a few Gyr ago 
in the present GMC-SC merger scenario. Recent numerical 
simulations have suggested that the MCs might have started 
their strong tidal interaction about 3 — 4 Gyr ago for a rea- 
sonable set of orbital parameters (Bekki et al. 2004b; Bekki 
& Chiba 2005) . Thus it is possible that the origin of SCs with 
DMSTOs may result from high GMC-SC collision/merging 
rates in the LMC due to strong tidal interactions between 
the MCs a few Gyr ago. 

Owing to continuous interaction between the MCs after 
the last dynamical coupling about 3 — 4 Gyr ago (e.g., Bekki 
& Chiba 2005), the LMC can retain an enhanced cloud-cloud 
(and thus probably GMC-SC) collision rate. This means 
that SCs with ages younger than 1.5 Gyr can also show DM- 
STOs, though strong observational evidence for the presence 
of DMSTOs in the CMDs of young SCs does not yet exist. 
Previous numerical simulations (e.g., Yoshizawa & Noguchi 
2003) have shown that the LMC and SMC can very strongly 
interact with each other about 1.5 Gyr ago - in this case it 
would be highly likely that a larger fraction of SCs with ages 
of 1.5±0.3 Gyr may show DMSTOs. This would be consis- 
tent with the recent observational results of Milone et al. 
(2008). In addition, the cloud-SC collision rate at the epoch 
of disk formation in the LMC might well be much higher 
than that at the present owing to the higher gas mass frac- 
tion and the higher degree of random motion at the time 
of disk formation. Therefore, it is very possible that the 
old GCs in the LMC may show DMSTOs in the GMC-SC 
merger scenario. 

It should be stressed, however, that recent proper mo- 
tion measurements of the MCs by the Advanced Camera for 
Surveys (ACS) on the Hubble Space Telescope (HST) have 
reported that the LMC and the SMC have significantly high 
Galactic tangential velocities (367 ± 18 km s _1 and 301 ± 52 
km s _1 , respectively), suggesting that the MCs may be un- 
bound from each other (e.g., Kallivayalil et al. 2006). In this 
case, LMC-SMC interactions about 1.5 and 3 — 4 Gyr ago 
are unlikely to have occurred, meaning that the formation 
of SCs with DMSTOs has little to do with any past inter- 
action between the MCs. Given that detailed modeling of 
the LMC-SMC orbits that includes the new proper motion 
data as well as factors previously ignored (e.g., a common 
halo and/or different circular velocity of the Galaxy) has not 
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yet been completed, it would be fair to say that at present 
the connection between the interaction histories of the MCs, 
and the formation of DMSTOs in SCs is not so clear. 

4.4 Relations to the Galactic GCs with 
abundance inhomogeneities? 

The origin of the observed star-to-star light element abun- 
dance inhomogeneities in the Galactic GCs (GGCs) has been 
extensively discussed both theoretically and observationally 
(e.g., Sneden et al. 1992; Norris & Da Costa 1995; Smith et 
al. 2005; see Gratton et al. 2004 for a recent review). The ob- 
served presence of star-to-star abundance inhomogeneities in 
less evolved stars on the main sequence and subgiant-branch 
(e.g., Cannon et al. 1998 for 47 Tuc) now strongly suggests 
that the origin of these abundance inhomogeneities is due 
to the early chemical evolution of GC-forming gas clouds. 
Stellar ejecta from AGB stars and massive stars has been 
considered to play a key role in the early chemical evolu- 
tion of GCs for the self-pollution scenario (e.g., Karakas et 
al. 2006; Prantzos & Charbonnel 2006). The early chemical 
pollution in GC-forming clouds results from stellar ejecta 
from massive AGB stars in the "AGB scenario" (e.g., Ven- 
tura & D'Antona 2005; Bekki et al. 2007) and from the stel- 
lar winds from massive stars in the "massive star scenario" 
(e.g., Prantzos & Charbonnel 2006). 

If the origin of the GGCs with abundance inhomo- 
geneities is essentially the same as that of the SCs with 
DMSTOs in the LMC, then age differences between the first 
and the second generations of stars in these GGCs may be 
as large as ~ 3 x 10 8 yr. This means that the massive star 
scenario can be ruled out, because the time lag between 
the first and the second generations is at most an order of 
~ 10 6 yr (Bekki & Chiba 2007; Decressin et al. 2007). This 
also means that gaseous ejecta from the first generation of 
AGB stars in a SC needs to be retained in the SC without 
star formation for ~ 3 x 10 8 yr and then converted into the 
second generation in the AGB scenario. Such delayed forma- 
tion of the second generation of stars implies that the second 
generation would form from the mixed gas of stellar ejecta 
from AGB stars with different masses. Thus, if only stellar 
ejecta from massive AGB stars (~ 6M ) can much better 
explain the observed C-N, O-Na, and Mg-Al anticorrelations 
(e.g., Ventura & D'Antona et al. 2008), the AGB scenario 
would have a serious problem in explaining the possible age 
different of ~ 3 x 10 8 yr between the first and the second 
generations. 

It would be equally possible that the origin of GCs with 
abundance inhomogeneities in the Galaxy is different from 
that of SCs with DMSTOs in the LMC, given that the phys- 
ical properties are different between these clusters in the two 
different galaxies (e.g., Mackey & Gilmore 2004). The im- 
plied high fraction of He-rich stars with Y > 0.3 in ui Cen 
and NGC 2808 (Piotto et al. 2005; Piotto et al. 2007) cannot 
be explained simply by the GMC-SC merger scenario: such 
He-rich stars would need to form from stellar ejecta either 
from fast-rotating massive stars (e.g., Prantzos & Charbon- 
nel 2006; Decressin et al. 2007) or from massive AGB stars 
(e.g., D'Antona et al. 2002). Possibly, the time-lags between 
formation of the first and the second generations of stars in 
clusters can be different between host galaxies. This might 
well cause a variety of differences in physical properties be- 



tween the first and second generations of stars in clusters 
belonging to separate host galaxies. 



5 CONCLUSIONS 

We have numerically investigated how GMCs evolve during 
GMC-SC merging in order to better understand the origin of 
the DMSTOs observed in intermediate-age SCs in the MCs. 
We have mainly investigated mass fractions of gas within 10 
pc of the remnants of GMC-SC mergers for variously differ- 
ent model parameters. We summarise our principle results 
as follows: 

(1) Gas initially in GMCs can be accumulated within 
SCs during GMC-SC merging so that high-density compact 
gaseous regions that are significantly flattened are formed 
in the central regions of the SCs. The mass fractions of gas 
(/ g ) within the central 10 pc of the merger remnants (i.e., 
the newly formed SCs) can be as large as 0.5 for models with 
s g ~ 1 (i.e., major mergers). 

(2) There is an optimum M s (or s g ) for a given set of 
model parameters, for which / g is a maximum. For example, 
the model with M B = M g (i.e., s g = 1) shows the maximum 
/ g in models with M g = 5 x 10 4 M Q and y g = 0.5i? g (=6 
pc). The remnants of major GMC-SC mergers are likely to 
show larger / g in the present study. 

(3) As long as GMC-SC merging occurs, / g does not 
depend strongly on y g and |u g |. The models with larger y g 
and |w g |, for which GMCs can interact with SCs without 
merging, show only very small f g (an order of 10~ 2 ). These 
results mean that merging is essential for the formation of 
new SCs with large f g . 

(4) If star formation is included in GMC-SC merging, 
compact, flattened star clusters composed of new stars can 
be formed in the central regions of the merger remnants. 
Since the simulated SCs-within-SCs appearances (i.e., dou- 
bly nested SCs) are not observed in the MCs, some later 
dynamical evolution processes need to transform the doubly 
nested SCs into normal-looking (well-mixed) objects with 
standard (King-type) density profiles. 

(5) The time scale of GMC-SC merging (t m ) in the LMC 
can be similar to the typical age difference of ~ 3 x 10 8 yr 
between the component stellar populations implied by ob- 
servations of DMSTOs in SCs in the LMC. However i m de- 
pends strongly on the number of GMCs, the sizes of GMCs, 
and the velocity dispersions of GMCs and stars in the LMC 
a few Gyr ago, all of which remain observationally unclear. 

Based on these results, we have pointed out that the 
observed possibly large fractions of the second generations 
of stars in SCs with DMSTOs in the LMC can be due to 
the past GMC-SC merging. We have also suggested that 
time lags between SC formation and the subsequent GMC- 
SC merging can be responsible for the possible age differ- 
ences between the first and the second generations of stars 
in the SCs with DMSTOs in the LMC. The adopted numer- 
ical code does not allow us to investigate whether the sim- 
ulated doubly-nested SCs can evolve into normal SCs with 
canonical radial density profiles due to internal dynamical 
processes. We plan to investigate this question in our future 
studies using the appropriate numerical codes (e.g., NB0DY4). 
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